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Abstract 

The Gauss-Kuzmin statistics for the triangle map (a type of mul¬ 
tidimensional continued fraction algorithm) are derived by examining 
the leading eigenfunction of the triangle map’s transfer operator. The 
technical difficulty is hnding the appropriate Banach space of func¬ 
tions. We also show that, by thinking of the triangle map’s transfer 
operator as acting on a one-dimensional family of Hilbert spaces, the 
transfer can be thought of as a family of nuclear operators of trace 
class zero. 


1 Introduction 


It is classical to think of continued fractions as iterations of the Gauss map 
on the unit interval and then to consider the underlying dynamical system. 
For example, the Gauss-Kuzmin statistics for continued fractions can be 
reduced to the study of the spectrum of the associated transfer operator 
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[la CHI Ea i3oi [31]. Most multidimensional continued fraction algorithms 
are iterations of a triangle A (though the Jacobi-Perron algorithm acts on a 
square). For background on many types of multidimensional continued frac¬ 
tion algorithms, see Schweiger’s Multidimensional Continued Fractions [51] 
and Karpenkov’s Geometry of Continued Fractions [2l]. Thus it is natural 
and standard to think of these as dynamical systems. This paper concen¬ 
trates on the transfer operator for the triangle map [91 [I3112I11291I3I] and 
then uses results on the spectrum of this operator to look at the triangle 
map’s Gauss-Kuzmin statistics. 

We will show that the transfer operator for the triangle map is 



Note that though this transfer operator looks similar to the transfer operator 
for the Gauss map: 



there is one significant difference, namely that in the triangle transfer oper¬ 
ator there is the term nx in the denominator, while the continued fraction 
transfer operator has the term n -|- x in the denominator. This difference 
is what prevents us from applying the standard bounds for the continued 
fraction case to the triangle case. It is this difference that to a large extent 
makes this paper not simply an easy generalization of earlier work. 

In this paper we find an appropriate Banach space of functions for which 
we can show that 1 is the largest eigenvalue, with one-dimensional eigenspace. 
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of the transfer operator. This can then be used to understand the Gauss- 
Kuzniin statistics for the triangle map. We then show, in analog to the Gauss 
map, that there is also a Hilbert space approach. We will see that in an 
appropriate sense the triangle map’s transfer operator can be thought of as a 
nuclear operator. Unlike the continued fraction case, we cannot be working in 
a Hilbert space of square-integrable functions on the two-dimensional domain 
A, as the natural eigenfunction with eigenvalue one is not square-integrable. 
Instead we will be thinking of the transfer operator as acting on a space of 
functions f{x,y) that are only square-integrable with respect to the variable 
y. Thus we will be acting on a family of Hilbert spaces parameterized by the 
variable x. 

There are two overall goals for this paper. First is to begin the spec¬ 
tral analysis for the transfer operator of the triangle map. Second, though, 
is to set-up the machinery for joint work with Ilya Amburg |3]. While the 
triangle map, to some extent, is just one example among many of possible 
multidimensional continued fraction algorithms, in mm it is shown that, 
by varying the triangle map in a natural way, a whole collection of both 
new and old multidimensional continued fraction algorithms can be gener¬ 
ated, called triangle partition maps. Further, in that work, it is shown that 
these triangle partition maps generate a family of multidimensional contin¬ 
ued fraction algorithms, called combination triangle partition maps, that in 
turn are shown to include many, if not most, known multidimensional con¬ 
tinued fractions. In [3], it will be shown that the transfer operators for half 
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of the triangle partition maps have behavior similar to that of the triangle 
map. This will allow us to apply the results of this paper to this other work. 
More interestingly, we will see that the transfer operators of the other half 
of the triangle partition maps have remarkably different properties. 

We would like to thank Ilya Amburg for a lot of help on this paper, in¬ 
cluding collaboration on the actual formulas for the Gauss-Kuzmin statistics. 
We would also like to thank L. Pedersen for many useful comments. 


2 Transfer Operator for Triangle Map 


Subdivide A = {(x,?/) G : 1 > x > ?/ > 0} into subtriangles 

Afc = {(x, |/)gA:1 — X — /c|/>0>l — X — (fc-l- 1)?/} 


Definition 1. The triangle map T : A —)■ A is defined by setting T = Tk for 

(x, y) G Afc, and in turn setting 

rj. , X [y l-x-ky\ 

\x X J 

The point (x, y) E A has triangle seguence (oq, Oi, 02, ...) if 

T^{x,y) G A„^. 


For more on the triangle map, see [91 na EH |29l [31] 

To dehne the triangle map’s transfer operator, we hrst compute the Ja¬ 
cobian: 


J(x, y) = det 


a /yN a /yN 

OX \xJ oy \x/ 


d /l—x—ky\ d /l—x—ky\ 
dx V X ) dy \ X ) 



1 



X 


1 


X 


3 ’ 






5 


for {x,y) e Afc. 

Next we need to find the inverses tk ■ —)• A of the triangle map, one 

for each non-negative integer k. These are 

tkix,y)= ( . ’ 1 I r I 

+ kx + y 1 + kx + y 

which can be checked by direct calculation,. 

By definition, for any differentiable map T, the transfer operator is 

This leads to 


Proposition 2. The transfer operator for the triangle map is 

^f){x,y) = J^TT—-( lit. I ’ 1 ^ ^ 

^ [1 + nx + yy^ \1 + kx + y 1 + kx + y 

This map is also be called the Perron-Frobenius operator or the Ruelle- 
Perron-Frobenius operator with respect to the Lebesgue measure dxdy, hav¬ 
ing the property that for all / G L^{dxdy), and all measurable subsets 
^ C lA^Tif)dxdy = fdxdy. For more on transfer maps in general, 

see Baladi [5] and for continued fractions, see losifescu and Kraaikamp [T8] . 

By direct calculation, for the function 



f(x,y) 


1 

x(l + y) 


we have that 


^f)(x,y) = f(x,y). 
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What we now must do is find appropriate vector spaces of functions on which 
the operator £ acts so that f{x, y) is the leading eigenfunction, meaning that 
we would know enough about the spectrum of C to be able to state that the 
largest eigenvalue is one, and that it has a one-dimensional eigenspace. 

3 Banach Space Approach 

We want to identify a Banach space of real-valued functions on A on which 
the transfer operator has leading eigenfunction with eigenvalue one. 

Let C(A) be the vector space of real-valued continuous functions on the 
triangle A. Note that since A is open, the elements of C(A) need not be 
bounded. 

Set 

B = {/ G C(A) : 3A G M such that \xf{x,y)\ < C,V{x,y) G A}. 

This is a Banach space under the norm 

11/11= sup \xf{x,y)\ 

(x,yEA 

Theorem 3. The transfer operator is a continuous linear map from V to 
itself. 

Proof. Linearity is immediate. We will show for all / G B that 

\\Cf{x,y)\\ < 3||/||. 

For any / G B, let A = ||/||, which means that for all {x,y) G A we have 

\xf{x,y)\ < C. 
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Then 

\xCf{x,y)\ = 

< 


X 


< 


E 

n={ 
oo 

E 

n=0 
oo 

E 

n=0 
oo 

^ \^(l + nx + yY 

OO 


V (1 

+ 

nx 

+ yf 

(_ 


X 


1(1 

+ 

nx 

+ yY 



X 


1(1 

+ 

nx 

+ yY 


Cx 




1 + nx + y^ 1 + nx + y 


f 


X 


n=0 


< Cx 

< Cx 

< 3C, 


{1 + nx + yY 
1 


0 V((l + I/)+ 


+ 1 


dt + 


{i + yf 


x{l + y) 


giving us that the transform map is indeed a continuous map of the Banach 
space V to itself. 

□ 


Before proving that the largest eigenvalue of C is one, with multiplicity 
one, we need a few lemmas. 

We have 

Lemma 4. For two functions f,gEV, suppose f{x,y) < g{x,y). Then for 

C{f){x,y) < C{g){x,y) 


for all {x, y) G A. 



Proof. This follows from looking at 


0 < g{x,y) - f{x,y) 


and then showing that 


0 = £(0) < C{g - f)){x, y) = C{g){x, y) - C{f){x, y), 


which is clear. 


□ 


We will also need 

Lemma 5. For any / G V, there is a positive constant B such that for all 
{x, y) G A, we have 

-B ,, , B 

x{l + y) x{l + y) 

Proof. Let / G Id be any element in V. Then there is a i? so that 

I ri B 

M 

(The reason for the 1/2 will be clear in a moment.) Then we have 

II ^ B 
^ 2x ^ x{l + y) ’ 

as desired. 

□ 

Theorem 6. The largest eigenvalue of C : V ^ V is one, with multiplicity 


one. 



9 


Proof. (The overall structure of this proof is similar to the corresponding 
result for the Gauss map. Again, what is new is hnding the correct Banach 
space and using the fairly recent result of Messaoudi, Nogueira, and Schweiger 
[29] that the triangle map is ergodic.) 

Let / G G be an eigenfunction of C with eigenvalue A. 

We know there is a constant B such that for all (x, y) E A 

_ Q Q 

-7:r -.—r < f{x,y) < —r- 

x{l + y) x{l + y) 

Then for all positive n, 


and hence 


Thus we need 


-B 

x{l + y) 


<C^^^f{x,y)< 


B 

x{l + y) 


-B 

x{l + y) 


<\V{x,y)< 


B 

x{l + y)' 


|A|<1. 


We now have to show that this eigenvalue has multiplicity one. To prove 
this from scratch would be hard, but using that the triangle map has been 
proven to be ergodic by Messaoudi, Nogueira, and Schweiger [29], the result 
follows immediately from Theorem 4.2.2 in Lasota and Mackey’s Probabilistic 
Properties of Deterministic Systems [23] which states 

Theorem. Let (X, A, y) be a measure space, S : X ^ X a non-singular 
transformation, and P the Frobenius-Perron operator associated to S. If S 
is ergodic, then there is at most one stationary density /* of P. Further, if 
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there is a unique stationary density /* of P and f*{x) > 0, a.e., then S is 
ergodic. 

Here P is our C. This gives us our result. 

□ 

It would be interesting to give an argument for this along the lines of the 
hrst few chapters in [23] . 

4 The Gauss-Kuzmin Distribution for Trian¬ 
gle Maps 

(This section is joint work with Ilya Amburg, and is a special case of chapter 
16 in | 2 ]). 

We would like to know the statistics behind a number’s triangle sequence. 
Knowing now that the transfer operator has leading eigenvector l/x(l + y), 
coupled with the work of Messoundi, Noguira and Schweiger [22] showing 
that the triangle map is ergodic, allows us to apply standard theorems to 
explicitly determine these statistics. 

Let {x, y) E A have triangle sequence (oi, 02 , 03 ,...). We dehne 

'■ cii = k and 1 <i <n} 

Pn,k{^,y) = -, 

n 

or, in other words, the percentage of the a* that are equal to k in the hrst n 
terms of the triangle sequence. Provided the limit exists, we set 

Pk{x,y) = lim pn,k{x,y). 

n—)-oo 
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We know that d/x = dxdy is an invariant measnre with respect to 

the triangle map T : A —)• A. (The extra 12 / 7 r^ is just to make the measure 
of the domain A be one.) Then in direct analog to Theorem 9.14 in izu , we 
have 


Theorem 7. For almost all {a, (3) G A 
Pk{a,f3) = [ d/r = 


12 


JA, j^^n^x{l + y) 

The proof is exactly analogous to that in 


dxdy. 


Theorem 8 . For the triangle map T, we have that 

6 Li 2 (i) + 121 og"( 2 ) 


F(0) = 1 - 


TT^ 


and, for k > 0, 
6 

P{k) = - 

TT^ 


Lio 




Lio 


{k + 2y 


,2/^ + 2 


+41og^(fc + l) -2log I 
Here Li 2 is the dilogarithm function. 


- 2 log{k{k + 2 )) log{k + 1 ) 


Proof. These are calculations, for which we used Mathematica. We have 
from the previous theorem that Pk{a,/3) is independent of {a, (3) G A for 
almost all elements in the domain. Then we have 


F(0) 


12 


/ a . Tr‘^x{l + y) 


dxdy 


12 


1 - 


7r2a:(l + y)^^ 

6Li2(i)+121og2(2) 


TT^ 


dx 
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and, for /c > 0, 
P{k) = 


12 


A, TT^x{l + y) 


1 — x 
k 


dxdy 


12 


_ ,.'13 A(l + j/) 


d?/ ) da: + 


1 

fc+1 


12 


_ ,.i5fA(l + s) 


dy ) da: 


TT^ 


Lio 


(k + iy 


— Lio 


^2/7 I 1 \ 01_2/^ + 2 


+4 log (A: + 1) — 2 log 


A: + 1 


(A:+ 2)2 

- 2 log(A:(A: + 2)) log(A: + 1) 


□ 


5 Attempts at nuclearity/ On various Hilbert 
spaces of functions 


For the Gauss map, Mayer and Roepstroff 123 showed that the transfer 
operator is nuclear of trace class zero. This section shows that nontrivial 
analogs hold for the triangle map. The difficulty in part stems from that the 
function l/a:(l + y) is not in the Hilbert space L‘^{dxdy), in contrast to the 
leading eigenfunction 1/(1 + x) of the Gauss map’s transfer operator being 
in L^(dx). Instead, we will think of l/x(l + y) as in L^{dy), treating the 
variable x as a parameter. 

On the positive reals, set 

dm{t) = — -dt. 

e* — 1 

Then set 

Vk{s) = 


(A: + l)! 












ek{t) 

Ekix,y) 
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= Llit) 

/CT(0(x,t)) = [ ^ 0(x,s)dm(5) 

Jo Vst e* - 1 

I _ 

0 (xy) = — e ^^(j)(x, s)dm(s) 

X Js=0 

poo 

(a(s),/S(s)) = / a{s)(3{s)dm{s). 

Jo 

(Here Ll{t) denotes the first Laguerre polynomial and Ji denotes the Bessel 
fnnction of order one.) 

Our goal is to show that 


ETf{x,y) = CT${x,y)) 

^ '^)Ct{4>{x, t))dm{t) 

= ^{<l>{x,s),yk{s))Ek{x,y) 

k=0 

After making a signihcant change of variables, we will see that our argu¬ 
ment mirrors that given in [27] , 

We start with functions f{x, y) with domain A for which there is a (j){x, s) 
that is in L^(dm(s)) such that f = (p cind such that Cxif) exists. Note that 
for 0(x, y) = 

^ = Hx,y))-, 
x{l + y) 

thus such functions exist and include the function l/a;(l -f y). 

We have 



ET{f){x,y) 




(1 + nx + yY 

T1=0 ^ ^ 


/ 


X 


1 + nx + y'' 1 + nx + y 
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E 


X 


{1 + nx + nY \1 + nx + y' 1 + nx + y 


Set 


w = 


1 + y 


X 


It is this change of variables that will allow us to use m 
Then we have 

1 1 


{1 + nx-\-yY x^{n + wY 

1 a; \ / 1 1 


1 + nx + y^ 1 + nx + y 
Then CT,^,if){x,y) is 

OO 

E 


x{n + w)^ n + w 


^x^(n + wY \x(n + w)^n + w 

n=0 ^ / \ \ / 


which is 


which is 


E 


/ + ■w^)) / e"+™0(a;, s)dm(s) 

;^^xYn + wY Jo 


7 —^-^e^0(x,s)dm(s). 

^^Jo ^oin + wY 


We want to rewrite 


E 


^ (n + w) 

n=0 ^ ' 


- 6 '^“ 1 "'^ 


to eliminate the dependence on n. 
We have 


E 


^ (n + wY 

n=0 ^ ' 


g n+iD 




OO -| OO 

V—^—V, 

^ in + wY ^ \ k\(n + wY 

n=0 '' ’ k=0 V ^ ^ 
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k=0 


k\ j in + 

n=0 ' ' 


From standard arguments on the Lerch Zeta function, we know that 


E 


r-oo 


^(n + M;)^+2 (/c + i)!7q i_e 


— t 




{k + 1)! Jq 




Thus 

OO 

E 


—^ (n + wy 

n=0 ^ ^ 


g n+it) 


“ "-s)* 


E 


k=0 
r*oo 




A:! J{k + l)\Jo 




\/st 


where Ji is the Bessel function of order one. 
Then 

^oo poo 


CT,t,f{x,y) = 




x^ 



'0 Jo 

r‘00 POO 



0 Jo 


Ji (2\^) ^ s)dm{t)dm{s) 

yst 

Ji(2\^) ^_^(' +y s)dm(t)dm(s) 

yst 


This is why we set 

lCT,t,{(t){x,t)) = 

We have 

k^T,Jix,y) = 


Ji(2\/st) 

\/it 


0(a;, s)dm(s). 



lm(t) 
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We want an understanding of the operator 
From m, we know that 


where L\{s) is a first Laguerre polynomial. 

We set 

efc(s) = Ll{s),r]k{s) = 

Both of these are actually functions of only one variable. 
Then, also from [27], we have 


fc =0 


/Cr,/.(0(a:,t)) = s),7]k{s))ek{t) 

as follows from 


Ji(2\/sf) 


\fsi 

) 


0(x, s)dm(s) 

0(a;, s)dm(s) 


s^e ^ 


k=0 


n + 1)! 

/ ^efc(f)?7fc(s)0(a;, s)dm(s) 

fc =0 

CXD „QO 

'^ek{t) / ?7fc(s)0(a;,s)dm(s) 

fc=0 

CXD 


A:=0 


The functions ?7fc(s) and efc(s) have already been studied in [27|, where it 
was shown that 

Vk{s),ek{s) e L‘^{dm{s)) 
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fc =0 


and, with respect to this inner product, 

C30 

\ek\ ■ \r]k\ < OO. 

Then we have 

^ [^(0(a;,s),%(s))efc(t)] dm(t) 

^(0(a;,s),r/fe(s)) 


^T,t,f{x,y) = 


,k=0 

1 


) 


k=0 

OO 



J e *( -’^'')efc(t)dm(t) 


Y,mx, s), Tikis))Ek{x,y), 


k=0 


where we set 


Ekix,y) = 



e ^ ^)efc(t)dm(t). 


We want to show that Ekix, y) G L^(d|/), again treating the variable a; as a 
parameter. (The calculation below will also yield that Ekix,y) ^ L'^idxdy).) 
We have 

t 


Ekix,y) = 



, x^ 



-t 


' l — x + y 


)L^(t) 


e* — 1 


-dt 



°° i-('^~^+y\ , -1 , , te * 

e"T—t --dt 


■X\ 

\ roo 





1 — e * 

OO 




—mt 


m=0 


f*oo 







m=0 


where s(m) = 
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1 \ ^ 

-.)T. ‘iiw- 

^ / rr,=n >^0 




m=0 

oo 


1 \ (n + l)(s(m) - 1)’' 
^2 j 


^' sim) 

m=0 ^ 


n+2 


by 7.414.8 in 


= (n + 1) 


1 ^ ^((i+^)_l)„ 


2 I / j / l+i/+ma;*^+2 

m=0 V X ) 


= (n 


i)E 


{1 + y + {m — l)xy 
{1 + y + mxY^'^ 


m=0 


As the above series grows like Yh there is a constant C independent 
of n snch that 


[n 




6 Conclusion 


Of course traditional continued fractions have many rich properties. It is 
natural to ask which of these properties have analogs for a given multidi¬ 
mensional continued fraction algorithm. This paper is the start of hnding 
the transfer operator analogs for the triangle map. In |3], similar analogs will 
be developed for triangle partition maps, which as mentioned, while being 
built in a natural way out of the triangle map, include many if not most 
multidimensional continued fraction algorithms. 

As mentioned in the introduction, transfer operator have been used for 
decades in the study of continued fractions, though people did not initially use 
the framework or rhetoric of functional analysis. The framework was made 
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explicit in the pioneering work of Mayer and Roepstroff IZHEH] and of Mayer 
[25l |26] , work that continues even today [36l EH [HI IH [20l [121 [HI El H EZl 
ElE]- For each of these papers, there are corresponding natural questions for 
multidimensional continued fractions. The corresponding questions should 
be non-trivial but interesting. 
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